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We investigate the dynamics of kinetically constrained models of glass formers by analysing the 
statistics of trajectories of the dynamics, or histories, using large deviation function methods. We 
show that, in general, these models exhibit a first-order dynamical transition between active and 
inactive dynamical phases. We argue that the dynamical heterogeneities displayed by these systems 
. are a manifestation of dynamical first-order phase coexistence. In particular, we calculate dynamical 

^ ' large deviation functions, both analytically and numerically, for the Fredrickson-Andersen model, 

D ^ the East model, and constrained lattice gas models. We also show how large deviation functions 

^ . can be obtained from a Landau-like theory for dynamical fiuctuations. We discuss possibilities for 

similar dynamical phase-coexistence behaviour in other systems with heterogeneous dynamics. 
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I. INTRODUCTION 



In this paper we describe in detail a theoretical method for the study of the dynamics of glassy systems [H, 0, S, 0] ■ 
T This approach is in essence a statistical mechanics of the trajectories of the dynamics, or histories, as it is based on the 
' O study of large deviation functions [sj — which can be thought of as generalized free-energies — of dynamic observables. 
^ In particular, we use the tools of Ruelle's thermodynamic formalism 0|i as applied to continuous time Markov 
chains [1], to study kinetically constrained models (KCMs) of glass formers @. In a recent letter ^] we showed 
using these methods that the dynamics of KCMs takes place on a first-order coexistence line between active and 
inactive dynamical phases, in accordance with previous suggestions flU. Here we expand significantly on Ref. ^10|, 
demonstrating in detail the existence of the first-order dynamical phase transition, and discussing the Landau-like 
approach [l^l that we use to characterise the dynamical phases, and the transition between them. The dynamical 
transition we find in KCMs is related neither to a thermodynamic transition, nor to a finite temperature (or finite 
psj ,. density) dynamical singularity. Our results, therefore, point towards a persp ective (llj on glasses which is distinct 
' from other approaches, such as the random first-order transition theory [Tsl. Il4l. [Tsl. Il6l|. frustration- limited domains 
'. [13 ) or mode-coupling theory [l^ . 

• The paper is organized as follows: in Section |TT] wc introduce our dynamical tools and the ensemble of histories 
OO ' ill which the dynamical phase transition takes place. In Section IIIII we describe the models that we will consider. 

■ We show the existence of a dynamical phase transition in Section IIVI comparing different models and establishing 
^ ' minimal conditions that are sufficient to ensure a dynamical transition. In Section |V] we discuss the ensemble of 
histories in detail, considering statistical properties of the active and inactive phases, and a dynamical analogue of 
phase separation. We summarise our results in Section IVTl and consider some open questions. 



II. DYNAMICAL TOOLS: THE s-ENSEMBLE 



A. Motivations 



In this article, we are concerned with fluctuations in dynamical observables such as the amount of dynamical 
activity in a glassy system, integrated over a long time t and over a large (but finite) system. To investigate these 
fluctuations, we consider statistical properties of the histories followed by the system. Ensembles of histories are 
central to the thermodynamic formalism developed by Ruelle and coworkers 0] (see for a comprehensive review). 
While thermodynamics is concerned with probability distributions over configurations of a large system, we will apply 
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the thermodynamic formalism to probability distributions over histories. We begin by discussing the physical content 
of the observables that we will consider. 

In the Boltzmann-Gibbs theory, the macroscopic features of large systems are characterised by determining the 
statistical properties (the mean value and fluctuations) of extensive observables, such as the energy or the number of 
particles. In a microcanonical approach, one considers the properties of a system with fixed total energy E. They are 
obtained from the counting factor 



n{E,N) = 



number of configurations 

. , (1) 
with energy E 



where N represents the size (the volume) of the system. In the large size limit {N — > oo), we define the entropy 
density s(e) = limjv^oo jf ln.Q{eN, N), which represents the relative weight of configurations with energy density e. 

In a dynamical context, we consider histories of the system between an initial time t and a final time t = t. 
Instead of considering the statistics of the energy i5, we will consider an observable A, that is extensive in the 
observation time t. The dynamical analog of 51 (i?, N) is the probability distribution of this observable 



fldy„(^,<) = 



fraction of histories with a given value of the 

(2) 

time-extensive observable A 



On a mathematical level, the choice of the observable A is somewhat arbitrary, although application of the thermo- 
dynamic formalism requires that the quantity j log f2dyn(i3.^i t) should have a finite limit for large times t. Subject to 
this constraint, the choice of the order parameter A is informed by physical insight: we should use an observable that 
reveals the essential physical processes at work in the system. For example, in non-equilibrium systems in contact 
with two reservoirs of particles, we might define A as the total particle current: the number of particles transferred 
from one reservoir to the other between times and t (see, for example, Refs [l^ [13]). In the context of glassy 
phenomena, we consider observables that measure the "activity" or the "complexity" of the history [1, [l^, [llj . 

Returning to the Boltzmann-Gibbs approach, it is useful to define the canonical ensemble through the partition 
function 

Z{(3,N) = Y,n{E,N)e-f"' (3) 

E 

which characterises a system at a given temperature /3~^. Within this framework, phase transitions can be identified 
from singularities in the intensive free energy, /(/3) = — limAr^oo In Z(/3, A^). The dynamical analog of this 
thermodynamic partition sum is 

ZA{s,t) = ^r!dy„(A,t)e-^^ (4) 

A 

where we introduced an intensive field s conjugate to A. This field will play a role similar to the inverse temperature 
f3. The dynamical partition function Zyi(s,t) is the central object of Ruelle's thermodynamic formalism. 

We have focused on the correspondence between the thermodynamic limit of large system size {N — > oo) and the 
long time limit {t — + oo) in Ruelle's formalism. In the following, we will consider systems for which the large time 
limit is to be taken at fixed system size: in some cases, we will then take a second limit of large system size N . If we 
consider systems with no thermodynamic phase transitions, then no singular behaviour arises on taking the limit of 
large N at fixed t. In this case, we expect the limits of large N and large t to commute, but this is is clearly not the 
case in general. 

B. Systems with Markov dynamics: statistics over histories 

1. Continuous time Markov evolution 

We now give more precise definitions of the quantities discussed so far, by reviewing the construction of the ensemble 
of histories for stochastic systems. We focus on continuous-time Markov dynamics (in this section, we follow @). The 
system is defined by a finite set of configurations {C}. Its dynamical evolution is defined by the rates W{C — > C) for 
transitions from configuration C to configuration C' . Thus, the probability P(C, t) of being in configuration C at time 
t evolves according to a master equation: 

dtP{C, t) = -r{C)P{C, t)+Y, W{C' ^ C)P(C', t) (5) 
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FIG. 1: A history of duration t is defined by a sequence of configurations Co ... Ck and a sequence of jump times 
ti, . . . ,tK- Between and tk+i, the system stays in configuration Ck. 



where 



ic) = J2w{c^c') 



(6) 



represents the rate of escape frora C. Equ.® is sufhciently general to describe kinetically constrained models, spin 
facilitated models, or lattice gases, with C representing the configuration of the whole lattice in each case. 

Starting from a configuration Cq at initial time t = 0, the system will experience a fluctuating number of changes 
of configuration ("jumps") between and t. We shall refer to the number of jumps as the "activity" and denote it 
by K. A history (or trajectory) consists of a sequence Cq . . . Ck of visited configurations, and a sequence of 
times ti, . . . ,tK at which the jumps occur (Fig.[T]). We stress that for a fixed observation time t, the number of jumps 
is a fluctuating quantity: it depends on the particular history followed by the system between and t. We refer to 
histories with many hops (large K) as 'active' histories and those with few hops (small K) as 'inactive'. 

We use the notation (O) for an average of the observable C, over histories of the system. We consider observables 
that depend on the entire history of the system, through the configurations visited and the time spent in each: that 
is, O — 0{Cq . . . Ck, ti . . . tK)- In general, we have 



K Cn-.-Cf 



K 



n ^(Cfe- 



.fe=l 



exp 



K 



^r(Cfe)(ifc+i - tk) 



k=l 



OiCo...CK,h...tK) (7) 



where the limits on the time integrals are ti > 0, tk < tk+i, and tK < tK+i — t; we use po(Co) to denote the 
probability distribution of the initial configuration Co- We use a compact notation for averages of this form: 



(O) = ^Prob [hist] O [hist]. 



(8) 



hist 



where Prob[hist] plays the role of a probability density in the space of histories. 



2. Time-extensive observables 



Having defined our system and its histories, we now turn to the choice of the time-extensive observable A. A simple 
choice of this observable will be the activity K. Each time the system changes configuration C C the activity K is 
incremented: K ^ K + 1. More generally, we can consider an observable A that is incremented at each jump, with 
the increment a(C,C') depending on the configurations before and after the jump. That is, for a given history with 
K changes of configurations 

K-l 

^[hist] - ^ a(Cfc,Cfc+i) (9) 

Again, we note that if a{C,C') = 1 then A is the activity K. 

To construct the dynamical partition sum, we start with a 'microcanonical' approach, classifying trajectories by 
their values of A. We generalise the probability P{C, t), defining P{C, A, t) as the probability of being in configuration 
C at time t, having measured a value A of the time-extensive variable between and t. Its evolution in time is given 
by the master equation 



dtPiC, A,t)^Y, ^(^' C)P(C', A - a{C', C),t) ~ r{C)P{C, A, t) 
c 



(10) 
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Thus, the probabihty of measuring a value A for the observable A in a history of length t is 

ndy,,iA,t)=Y,PiC,A,t) (11) 

c 

which we identify as the quantity introduced in ([2]). 

3. Canonical description: evolution in the s-ensemble 

We have defined the distribution ridyn(^: t) that is the analog of the microcanonical counting factor ^{E, N). We 
now introduce the analog for the canonical (Boltzmann-Gibbs) ensemble, parameterized by a field s. This involves a 
modification to the statistical weight of each history: 

Prob[hist] ^ Prob[hist]e-"'^['''"*l (12) 

Thus, in the 's-ensemble', averages of observables O are given by 

^ YJ^) 1^ ^5MProb[hist]e-^['^'-l = (13) 

where 

Za(s, t)^Y^ Prob[hist]e-"^[h'"'l = (e""^) (14) 

hist 

is the dynamical partition function, introduced in ((31) . (The subscript A of Za serves as a reminder that the field s 
is conjugate to A.) 

Averages in the ensemble with s = correspond to the steady state averages of O. A 'priori, this is the only 
physically accessible ensemble. Positive or negative values of s favor histories with non-typical values of A. For our 
purposes, working in the s-ensemble is simpler than considering ensembles with fixed values of A. We take the Laplace 
transform of P(C, A, t) with respect to A: 

PA(C,s,t) = ^e-^^P(C,A,t) (15) 

A 

From (fTU)) . the equation of motion for P4(C, s, i) is 

atP4(C,s,t) = ^e-^"('^'''^)M^(C'^C)PA(C',s,t) -r{C)pA{C,s,t), (16) 
c 

or, in an operator notation, dtpA — ^ aPa , where Wyi operates in the space of configurations {C}. Its matrix 
elements are 

(W^)c,C' = ^(^' ^ C)e-^"(c'''^' - r(C)^c,C'. (17) 

Some properties of the operator W are discussed in appendix[Xl Equ. (jA2[) states that -Pa(C, s, t) behaves in the large 
time limit as PA(C,s,t) ~ i?o(C, s)e*'^'*^'*^ where V'a(s) is the largest eigenvalue of Wa and i?o(C,s) is the associated 
eigenvector. Thus, for large times, 

ZA(s,t)=^P4(C,s,i)^e*'^^(^), (18) 
c 

and we will refer to V'a(s) as (the negative of) the dynamical free energy per unit time. Summing Eq. (|15p over C, 
probability conservation implies Z^(0,t) = 1, so that V'(O) = for all stochastic systems. 
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4- Large deviation functions 

In the Boltzmann-Gibbs theory, entropy and free energy are related through a Legendre transform (as can be seen 
from ^ or [2l|) which provides a hnk between microcanonical and canonical ensembles. We have already defined the 
function tpx{s), which is the dynamical analog of the free energy density /(/?). The dynamical analog of the entropy 
density s(e) is 

7r(a) = lim - log ildyn(a^, i) (19) 

which determines the large-i scaling of the probability of observing a value at for the observable A. 

For large times, the sum in Q is dominated by the maximum of ridyn(^, t), so that 7r(a) and i^A{s) are are related 
through a Legendre transform: 

i^Ais) = max (7r(a) — so) (20) 

If the function 7r(a) satisfies vr"(a) < 0, it can be obtained from the inverse transform 

7r(a) = min (?/'^(s) + sa) (21) 

Physically, the quantity 7r(a) describes the large fluctuations of A. It is maximal at the most probable value of a, 
which is the mean value oi A/t, in the limit of large time t. Gaussian fluctuations oi A/t are described by the quadratic 
approximation of 7r(a) around its maximum. Expanding 7r(a) beyond quadratic order gives information about non- 
Gaussian fluctuations of A/t^ which are referred to as large deviations 0. Alternatively, one may characterise 
these fluctuations through ij)A{s), since the cumulants of A are obtained from the derivatives of i^A{s) through 

limf^oo \ {AP)c — (—1)^ dsp^^"* ' where, as usual, {Ap)c is the p-th cumulant of A. 



5. Time-extensive ohservahles varying continuously in time 



In addition to time-extensive order parameters of the form given in ([9]), we also consider those of the form 

B[hist] = Y,{tk+i - tk)h{CK) = / dt' b{C{t')), (22) 
fc=0 -^0 

where we introduced a configuration-dependent observable 6(C). In the sum over fc, we define to = and tK+i = t so 
that the time spent in configuration Ck is simply tk+i — tk- In the integral representation, we have represented the 
trajectory by a function C{t') which takes the value Ck for tk < t' < tk+i- The time- integrated energy of the system 
is an observable of the form B, in which case b(C) is simply the energy of configuration C. Then, defining P{C,B,t) 
by analogy with P{C,A,t), we have 

9tP(C,i?,i)=X]W^(C'->C)P(C',i?,i) -r{C)P{C,B,t)-b{C)-^P{C,B,t) (23) 

We define an s-ensemble associated with the observable B through 

Prob[hist] ^ Prob[hist]e-"'^[h'"'l (24) 

Then, repeating the analysis of Section fll B 3[ the analog of iPa{s) is (t>B{s) = lirtit^oo jln(e^*-^). This quantity is 
equal to the maximal eigenvalue of an operator Wb , whose elements are 

(Wb)c,c' = ^ C) - [r{C) + sb{C)] Scfi' (25) 

In the following, we concentrate our study on time-extensive variables of type A. Some connections between s- 
ensembles parameterized by observables of types A and B discussed in Appendix [Bl 
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6. Variational approach for %1)a (s) 

The models considered in this work have dynamics which obey detailed balance with respect to an equilibrium 
distribution Peq(C): that is, Pcq{C)W{C C) = Peq(C')VF(C' — > C). This allows us to derive a variational bound on 
the dynamical free energy tpxis)- To achieve this, we symmetrise the evolution operator W^, defining through 
the similarity transformation {Wk)c,C' = Pcq'^^ iC)(W k)c,C' Pcq^ (C) ■ Hence, 

(Wk)c,c - [W{C ^ C')W{C' -^C)]^ - r(C)<5cc' - (Wa-)cc' (26) 

Since and W^^ are related by a similarity transformation, their eigenspectra are identical. We therefore use a 
variational principle (valid for any symmetric operator) to determine their common maximal eigenvalue: 

, , , Ecc ViC)iWK)c,C'V{n {V\Wk\V) 

^"^^^ ^ {^fc^i eTWP ^ T.T ^ywr ^''^ 

At s = 0, the maximum is achieved for V{C) — Pcq(C)^^^, and V'(O) = 0, as required. 

Interestingly, the quantity to be maximised in (j27p has a physical interpretation. For any history of the system, the 
fraction of time spent in each configuration C defines a quantity known as the experimental measure. As we discuss 
in appendix [Cl Donsker-Varadhan theory relates the probability of observing a particular experimental measure to 
an expectation value of the form (T^|Wi^|y). In Section|Vl we will use these results to investigate fluctuations in the 
s-ensemble. 



III. MODELS AND ORDER PARAMETERS 



A. Kinetically constrained models: FA, East, TLG and KA models 

Kinetically constrained models 0, il, [13, il 111, [13, [11, |M H ill, 111 [11, il] are simple lattice models of 
glasses which can account for a large range of dynamical phenomena associated to the glass transition. This includes: 
super- Arrhenius temperature dependence of timescales, non-exponential relaxation, spatially heterogeneous dynamics, 
transport decoupling, and aging and memory effects. The thermodynamic properties of KCMs are very simple, and 
their non-trivial features arise from dynamical rules which forbid or favor some transitions, while maintaining detailed 
balance with respect to a trivial equilibrium distribution over configurations. For a review on KCMs see 0]. 

We first consider models with binary spins rii = 0, 1 where i = 1, . . . , N are the sites of a lattice. In spin-facilitated 
models, sites with rii — 1 represent excitations, which promote local activity. The models evolve by single spin-flips, 
which occur with rates 

e/3(".-i) 

Win, ^ 1 - n.) = C\{{n,})j-^ (28) 

where P is the inverse temperature, and the kinetic constraint enters through the function Ci{{nj}), which is a function 
of the neighbors rij of i, but does not itself depend on Ui. In this case, it is simple to verify that the model obeys 
detailed balance with respect to the equilibrium distribution 

g-/3ni 

Peq{{n^}) = n 

i 

and that the excitation density is c = {rii) = (1 + e^)^^- 

In the one-spin facilitated Fredrickson- Andersen (FA) model [121 , Ci — 1 if any of the nearest neighbors j of 
i are in the excited state, Uj = 1; otherwise Ci — 0. We also consider the three-dimensional variant of the East 
model [H,[13,[3H| in which Ci = 1 for site i = (x, y, z) if at least one of the sites (x— 1, y, z), {x,y— 1, z) or (x, y,z—l) 
is in the excited state; otherwise d = 0. 

In addition, we consider lattice gas models [^,[13, [IE], in which particles move from site to site, with at most one 
particle per site. Sites which are occupied have rii = 1, and unoccupied sites have rii = 0. Particles move between 
sites i and j with rate Cij{{nk}) so that the model has a conserved density p = J^i^i- The rate Cij{{nk}) 
is non-zero only for nearest neighbor sites i and j, and it is independent of rii and rij. Thus, equilibrium state has 
a trivial distribution: all configurations with density p have equal probability. As an example of such a model, we 
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consider the two- vacancy facilitated triangular lattice gas, or (2)-TLG [Ss*], which is defined on a triangular lattice, 
with a constraint which is equal to unity if the two common nearest neighbors of sites i and j are vacant, and 
zero otherwise. Similarly, the (2,2) variant of the Kob- Andersen (KA) lattice model [13] is defined on a square lattice, 
with Cij = 1 if at least one neighbour k ^ j of site i has Uk — Q and at least one neighbour k' ^ i oi site j has nk' — 0. 
Otherwise Ca = 0. 



B. Reducibility of KCMs and sums over histories 

The construction of the s-ensemble in Section [IT] assumed that the system of interest has a single steady state to 
which it converges in the long-time limit. For finite-sized stochastic systems, this convergence is ensured as long as 
the dynamics are irreducible that is, it must be possible for every configuration of the system to be reached from 
every other configuration. For KCMs, this is not the case in general. For example, in the FA model, there are no 
transitions either into or out of the configuration with no excited sites {rii = for all i). For the models considered 
here, these states are usually considered to be irrelevant because they have a contribution to the Gibbs measure that 
vanishes exponentially in the thermodynamic limit, at all temperatures T > 0. 

However, when considering large deviations, these states may become relevant. In order to ensure convergence to 
a single steady state, we define our unbiased measure over histories as in ([7]), with a distribution of initial conditions 
po(Co) that is non-zero only for configurations in the largest irreducible partition of the dynamics. That is, we do not 
allow the system to occupy configurations that cannot be reached from representative configurations taken from the 
relevant Gibbs ensemble. For the FA model and East models, this simply means that the system may not occupy the 
configuration which has no excited sites. Practically, this means (for example) that the maximum in the variational 
expression (j27p should be taken with the constraint that V(C) is finite only for configurations in the largest irreducible 
partition. 

Instead of restricting initial conditions in this way, one could instead consider large deviations in a subsystem of 
size N that is embedded in a larger system of size N' ^ N: this was the approach taken in [llj. As usual, we expect 
these approaches to be equivalent in the limit of large system size N. 



C. Kinetically constrained models : bosonic and mean- field variants 

It is convenient to define a 'bosonic' variant of the one-spin facilitated FA model [1^, H^l, in which may be any 
integer greater than or equal to zero. We take 

Win, ^ n, + f) = C,({n,})e-^, Win, ^ n, - I) = ai{n,})n, (30) 

where Ci is again independent of 7ii so that 



ni . 



where c = (rii) — . For the bosonic FA model in finite dimension, we take Cii{nj}) =^ "-j: where the sum is 
over the nearest neighbours j of site i. 

For this bosonic model, it is convenient to use the Doi-Peliti representation [38| . We define bosonic operators a^, 
and hi — alui, with [a^, Oj] = 6ij and a vacuum state |0) for which ai\0) — for all i. The Doi-Peliti representation of 

the operator W is defined by (W)cc' — (0| Yli TTT ^ni('^i)"'|0) where the configurations C and C have occupations 
{ni} and {n'} respectively. In the s-ensemble, we are interested in the operator Wa defined in (|17p . In the case where 
the observable A is the activity K, we have 



l^e ''ical + ai)hj + e '*(caj + aj)hi — 2ihifij + c) (32) 



where the sum runs over (distinct) pairs of nearest neighbours. 

In addition, it is useful to consider a mean-field variant of the FA model, in which the facilitation function of site 
i depends symmetrically on the state of all sites. That is, Wini ^ ni + 1) = N^^^-nje'~^^ and Wini + 1 ^ 
Ui) = N^^ TijTii, which satisfy detailed balance with respect to ((3T|) . In the Doi-Peliti representation, the master 
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operator is simply 



,(FA,mf) 
'if 



(2iV) ^ ^ l^e *aj(ca] + aj)a^ + e *a](caj + ai)aj - 2(aja]ajai + c) 



(33) 



Due to the symmetry between sites, the properties of the model can be obtained from a single co-ordinate: the total 
number of excitations ntot = X)i whose equilibrium distribution is Poissonian with mean cN. In this co-ordinate, 
the master-like equation (|16p has a closed form, and the matrix elements of the relevant operator are simply 

(y^K^''^^^)ni^„n,^, = cntot(e~"5,H„,+i,„;^^ - '5nt„t,nu) + ^('^tot " 1) (e^"(5„^„^_i,„j_^^ - <5„,„,,„;^J (34) 



D. The A-model and the AA model 



It will be useful to compare the FA model with two other models, which we call the A and AA models. These 
names are motivated by the schematic representations of their fundamental processes, as A ^ and A + A ^ (d. Here 
we have used an alternative notation to avoid confusion with the observable A used to define the s-ensemble. 

We define the A-model and its bosonic variant by removing the kinetic constraints from the FA model: that is, 
Ci{{nj}) = 1, independent of the state of the system. In this model, excitations are created and destroyed singly, 
independent of site. The A-model has the same equilibrium distribution as the FA model, but its large deviations can 
be solved exactly. 

We also compare the FA model with a model in which particles appear and annihilate (AA) in pairs [s^. This 
so-called AA model is related to a variant of the FA model, through a similarity transformation that connects their 
master operators [37j . The AA model is defined for binary spins = 0, 1. In this model, the excitations move between 
adjacent sites with rate D, and appear and annihilate in pairs with rates k and fc' respectively. Schematically, we 
write 

0,l, <^l,Oj, 1,1, ^0,0j, O.Oj^ia, (35) 

for neighbouring sites i and j. The equilibrium state of this model is of the form ((29|) . with = ^/ k' /k. It is also 
convenient to consider bosonic and mean-field variants of this model, defined analogously to their FA counterparts. 
In the bosonic AA model, we generalise to Ui > 0, using rates 

{ni,nj) {rii - l,nj + 1), (n,;, n,) (n^ - 1, 71, - 1), (n^, n,) ^ (n^ + 1, + 1) (36) 
which obey detailed balance with respect to ([3T|) with k' = ke~'^^ as before. In the Doi-Peliti formalism, we have 

W^"^^ = e^^[fc'a^a] -I- ka.a^ + D{a\aj + a]ai)] - [kh.hj + ck' + D{n, + fij)] (37) 
(u) 

In the mean field variant of the AA model, diffusion occurs between all pairs of sites [i ^ j), with rate (D/N); pair 
creation and annihilation processes occur with rates (k/N)ninj and {k' /N) for all pairs of sites i ^ j; and we also 
allow for on-site pair creation and annihilation: ^ ri^ ± 2 with rates (k' /N) and {k / N)ni(ni — 1). In the Doi-Peliti 
representation, the master operator is 

W^^'"" = 2^ E [e"'(fc'«I4 + ka^a,) - {ka\a]a,a, + ck')] + ^ E [e"'(«I«^- + " ("» + ^i)] (38) 

For finite systems, the restriction to i j in the diffusion term means that the master-like equation cannot be written 
in terms of the single co-ordinate rttot, except at s = 0. However, in the limit of large- A'^, this single co-ordinate is 
sufficient, and the master-like operator for this co-ordinate reduces to 

(W^^''"^Viot,«tot = k'N{zSn,,,+2,nl^^ - Sn,at,n[J + k^{ntot " 1) (z5„^^^_2,ni^, ~ Kat^^J + D{z - VjUtot (39) 

with z — e~'^ . 



E. Relevant observables 



We now discuss the observables that we will use to define the s-ensemble, and those that we will use to characterise 
trajectories within that ensemble. 
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1. The activity K and the complexity Q+ 

We have aheady defined the activity K, which counts the number of changes of configuration in a dynamical 
trajectory. In the context of dynamically heterogeneous systems such as glass-formers, the local activity can be used 
to distinguish mobile and immobile regions of the system. The large deviations of the extensive activity K are used 
to characterise trajectories which are more or less mobile than average. 

We note that K is of the form given in ^ with a{C',C) = 1, so the properties of the relevant s-ensemble are 
encoded in the operator 

(W/Occ = ^(^' ^ ^) - KC)5c,C' (40) 

Systems with dynamical heterogeneities are likely to present a wide distribution of very different histories. One 
way of characterizing this diversity is provided by the dynamical complexity of the histories pTI . |40| . l4l| . In the 
context of dynamical system theory, this quantity is called the Kolmogorov-Sinai entropy ^ . It provides one with the 
information content of the history and is defined as the logarithm of the probability of the history. As discussed in , 
the appropriate generalisation of this approach to systems with Markov dynamics is to consider the entropy associated 
with the measure over sequences of configurations Co . . . — > Ck @ • This amounts to performing a coarse-graining 
in time: it means that the information associated with the time intervals between changes of configuration is ignored 
when calculating the complexity. The definition of the dynamical complexity is 

which is of the form given in Thus, we define a dynamical partition sum 

Z+{s,t) = {e~^Q+). (42) 

The corresponding dynamical free energy is V'-f- (■?) — limt_»oo j In (s, t) which corresponds to the topological pressure 
of dynamical system theory. The analog of the Kolmogorov-Sinai entropy /iks is 

/iKS = - lim i(g+) = ^V+(s) (43) 
t-+oo t as 

which provides a measurement of the dynamical complexity of the histories in the steady state. In the examples 
of glass formers we will study below, the dynamical ensembles given by K and are qualitatively similar: we 
concentrate on the activity K for simplicity. 



2. Fluctuation theorem in the s-ensemble 



The Gallavotti-Cohen relation holds also in the quasi-stationary state at fixed value of K (or s), and therefore 
the fluctuation-dissipation theorem holds there as well. In order to see this, we parallel the reasoning presented by 
Lebowitz and Spohn in 42], and we construct the operator governing the dynamics not only at fixed value of the 
activity K but also at fixed value of the entropy current Qs — '^n=o !y(c"+^^'c'^) ' which, in terms of the variables 
s and A conjugate to the activity K and the entropy current Qs respectively, leads to the following pseudo-evolution 
operator, 

(w(s,X)) =e-'W(C' ^C)^~^W{C' ^C)^ ~r(C)5cc' (44) 

whose property W(s, A)''' = W(s, 1 — A) ensures that its largest eigenvalue ^ verifies ip{s, A) = ^p{s, 1 — A). For system 
with particle conservation, and subject to a field driving the system out of equilibrium, we note that the entropy 
current Qs is directly proportional to the total current of particles flowing through the system [i^. In that case, the 
generalized symmetry ip{s, A) = "(pis, 1 — A) implies a fluctuation-dissipation like relation in the s-ensemble. 



3. Order parameters within the s-ensemble 

As well as using the observables K and to define s-ensembles through , we also characterise the s-ensemble 
by using two other order parameters. For spin-facilitated models, we consider the mean excitation density: 

^-(^)-/i-^(/d-E"^W>. (45) 
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For lattice gas models, the particle density is specified by the initial conditions, so we require a different order 
parameter. The average activity is given by ■^{K)s. One can also consider the mean escape rate r{C) which depends 
only on configurations of the system. Again, we time-average this quantity along the trajectories, and divide by the 
system size N, defining 
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rK{s)^lim—[ I drr(C(r))) (46) 

IV. DYNAMICAL TRANSITIONS IN MODELS OF GLASS-FORMERS 
A. Existence of a transition in KCMs: variational bounds 

It is clear from their equilibrium distributions Pcq(C) that KCMs have no phase transitions at any finite temperature. 
That is, their thermodynamic free energies are analytic functions of temperature (or chemical potential). However, 
we now show that in the limit of large time t and large system size N, the dynamical free energy density N~^tpK{s) 
has a singularity at s = 0. To be precise, the dynamical free energy has a discontinuous first derivative with respect 
to s, so we interpret this singularity as a dynamical analog of a first-order phase transition. 

The proof of such a transition is based on the escape rates r(C) from the configurations of the model. We establish 
two bounds on ip{s). Firstly, the number of jumps K is non-negative, so Eq. (|14p implies that Zx{s,t) is a non- 
increasing function of ,s. Thus ipKis) is also non-increasing. Further, ■0-R'(O) = 0, so we have 

^Kis) < 0, ■s>0 (47) 

Secondly, we can use the variational result ([?/)) with V{Ci) — 1 for just one configuration Ci, and V{C) = otherwise 
to establish 

V'(s) > -mjn[r(C)] (48) 

for all s. For our purposes, the most important property of the kinetically constrained models defined above is that 
they have 

lim iV-^min r(C) = (49) 

N^OD C 

This can be established by explicit construction. In the FA and East models, we simply consider a configuration 
containing exactly one excitation, which has escape rate 2dc in the FA case and dc in the East model, where d is the 
spatial dimension (in the bosonic variants, these rates are 2de~^ and de~^). In the (2)-TLG, all of the particles in 
the model can be arranged in a single compact cluster, in which all but six of the particles are unable to move: this 
configuration has r(C) = 6. For the KA model, a similar construction leads to configurations with r(C) = 4. Thus, 
combining (|47ll49p . we have 

lim 1v.k(s) = 0, s>0 (50) 
Recalling that (K) = t{d/ds)ipx{s), We define the mean activity per site per unit time as 

IC{s) ^ hm 1^ (51) 
^ ' N^oc N ds 

and we can see that 

/C(s) =0, s > (52) 

Further, from Eq. (jB9p . we have K{0) = t{r) = tj^c Pcq{C)r{C). Since the distributions Pcq{C) have simple forms in 
kinetically constrained models, this quantity can be calculated explicitly: the limit /Ccq — limAr_,oo Pcq{C)r{C) 
is finite and positive for all the models that we consider. Finally, it follows from (|14p that /C(s) is non-increasing, so 
that 

/C(S) > /Coq, s < (53) 

with /Ceq finite. Eqs. ([5^ and ([55]) establish the discontinuity of /C(s) at s = 0: in the limit of large system size, the 
dynamical free energy has a discontinuous first derivative which we refer to as a first-order dynamical phase transition. 
We have established the existence of such a transition in the FA, East and (2)-TLG models, in all dimensions and for 
all finite temperatures [and for all finite densities p in the (2)-TLG]. That is, the simple phase diagram shown in fig [5] 
is generic to all of these models. 



11 



Active 



Inactive 



CriticaliTy 



0.6 
0.5 
0.4 
0.3 



CO 

-IS 0.2 



0.1 




. active 


inactive 









-0.2 



-0.1 





S 



0.1 



0.2 





0.06 




0.05 


CO 


04 




0.03 


is- 




J< 


0.02 


N 






0.01 









-0.01 



\ active 


inactive 







-0.2 



-0.1 





S 



0.1 



0.2 



FIG. 2: (Left): Generic 'dynamic phase diagram' for spin-facilitated KCMs such as the FA and East models. There is a 
dynamical phase coexistence boundary at s = 0, for all finite temperatures. The boundary ends in a dynamical critical point at 
s = T = 0. For the (2)-TLG and KA models, the picture is identical if the temperature T is replaced by the fraction of vacant 
sites, 1 — p. (Center): Variational estimates for the activity per site /C(s), in the bosonic FA model, at c = 0.25. For s > 0, the 
result lC(s) = is exact in all dimensions, from Eq. (|52|) . For s < 0, the dashed line shows the lower bound obtained from ([53}, 
while the solid line is the variational estimate /Cvar(s) — — ■^(d/ds)i/)var(s), obtained from (|58p . As discussed in the text, the 
solid line gives the exact result for the mean-field variant of the FA model. (Right) Again, we show the exact result %Pk{s) = 
for s > 0, together with the variational lower bounds (|58|) (solid line, exact for the mean-field variant) and (148 |l (dashed line). 



B. Variational free energy for the excitation density pk{s) 

The analysis given above establishes some minimal conditions that are sufficient for the existence of a first-order 
transition. For a more quantitative analysis, it is useful to use a specific variational distribution in (j27p . We consider 
a general bosonic KCM with single spin-flip dynamics, and we define a distribution of the excitation numbers rii that 
is independent of the site «, and parameterized by a mean density p: 



^p({-.})-nv^ (54) 

From (j27p . we therefore have ip{s) > —N mirip J-k{p, s) with 



^k{p,s) = N ^y-^ (55) 

The value of p which minimises J-k{p,s) is denoted by /Cvai(s)- It represents a variational estimate for the order 
parameter pk{s): if the variational bound (j27p is saturated then \V) is an eigenvector of the symmetrised operator 
Wk, and it follows from ((Xl4)) that pk{s) = Pvar(s)/(1 - e"''^). 

For the bosonic FA model, it is straightforward to calculate Tk{p, s). The only subtlety is that we must explicitly 
exclude the state with no excitations from the inner products, as discussed in Section FlII Bl In the Doi-Peliti formalism, 
our choice for V{{ni\) renders this calculation very simple: in terms of the symmetrised operator W^"^\ we have 

(0|e- VpS:. w^^e- VP S:. 4 \q)^-pN 
Tk{p,s)=N _ . (56) 



Hence, 



T,ip,s)^2 f^[z'^:;/^% . (57) 



Minimising over p, we find liniTv^oo N > ipv&iis) with 



■'/'var(s) = yPvar(s)[/3var(s) " c] (58) 



and 



0, s > 

(c/8) (9e-2^ - 4 + 3e-V9e-2'* - , 
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FIG. 3: (Top left) Variational free-energies for the FA model (constrained dynamics, solid line), and the A-model (unconstrained 
dynamics, dashed line) at s = 0. Both models have the same thermodynamic free energy. However, the dynamical function 
T(yp, s) reveals that the FA model has two dynamical phases while the A-model has only an active phase. (Top right) The 
variational estimate pvar(s), for the FA model (solid line) and the A-model (dashed). For the mean-field FA model, pvar(s) 
coincides with Pk{s) in the limit of large system size A*'; for the A-model, Pk{s) = Pvar(s) always. (Bottom) Dependence of 
the variational free energy on the field s, in the FA model. At the phase coexistence condition, s = 0, the free energy has 
degenerate minima. For finite s, either the inactive or active phase is preferred. 



(Within this approach, we obtain Pvar(O) — c by minimising J-{p, s) at fixed system size N, and then taking N — > oo.) 
The bound on ipxis) and the corresponding estimate of /C(s) are shown figO The variational estimate for Pk{s) and 
the variational free energy Tk{p, s) are shown in Fig[3l 

So far, we have used Eq. (|27p to obtain variational estimates for 'ipxis) and pk{s) for the FA model in finite 
dimension. For the mean-field variant of the FA model, it can be shown that these variational estimates are exact, 
in the limit of large system size N . (The factor 2d that appears in !F{p, s) is simply an arbitrary rescaling of time in 
the mean- field model. Our definition of the mean-field model requires that we set 2d = 1.) That is, the difference 
between the variational ansatz of (|54p and the dominant eigenvector of W k vanishes at large N . Mean-field models 
are discussed in more detail in Section fVl below. 

It is useful to compare these results for the FA model with the bosonic variant of the A-model, for which it can 
be easily verified that the large deviation function iPk{s) coincides with the variational bound V'var(s), even for finite 
system size N. In that case, we have 

T{p,s)^c + p-2e-'^ = c(e-2^ - 1), ^^^(s) = ce'^^ (60) 

Thus, while constrained FA model and the unconstrained A-model possess the same equilibrium distribution Poq(C), 
and hence the same static free energies, their dynamical free energies show dramatic differences. For large systems, 
the FA model exhibits a dynamical phase transition, while the A-model does not. See fig.[31 

The presence of the dynamical first-order transition in the FA model is intimately connected to the two minima in 
^k{Pi s). As shown in Fig. [31 these two minima represent an active phase, with p c and an inactive one with p ~ 0. 
The global minimum of !F{p^ s) lies in the active phase for s < 0, while it lies in the inactive phase for s > 0. For 
s = 0, one must consider carefully the limit of large system size N: we have Tk{c.-, 0) = while the inactive minimum 
occurs at pvar — 0{N^'^), where the value of the variational bound ^var is positive. Thus, the global minimum of 
J-k{p, s = 0) occurs at the active state density. However, any s > is sufficient to drive the system into the inactive 
phase. The effect arises because of two non-commuting limits: when minimising J-{p, s), taking the limit s — > before 
the limit of large N results in active behaviour; on the other hand, taking the limit of large N followed by a limit 
s ^ 0"*" leads to the inactive phase. 

We note that the dynamical phase transition in the FA model requires a limit of large system size {N oo) as 
well as a limit of long trajectories (t — > oo). To keep our methods well-defined, we excluded the configuration with 
no excitations from the initial conditions (recall Section fill B[) . We emphasise that we have proved the existence of 
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a dynamical phase transition in an irreducible model, with no absorbing states (this can be compared, for example, 
with phase transitions in the directed percolation universality class 14 3|). 

C. Numerical results 

1. Cloning method 

We now present some numerical computations of the dynamical free energy iPk{s) in KCMs. From this can be 
obtained from the large t limit of the equilibrium average (e~'*^). However, direct calculation of this average requires 
a computational effort that scales exponentially with t: the average is dominated by rare histories lying in the tails 
of the distribution of K. In dynamical systems [31 and in discrete time Markov processes [l^, this problem can be 
avoided by using a cloning method similar to that used in quantum mechanical Diffusion Monte Carlo algorithms [45j . 
This method was generalised to continuous time Markov processes by Tailleur and Lecomte [20]. Here, we briefly 
summarize the algorithm for obtaining dynamical free energies. 

The function iPa{s) is obtained as the largest eigenvalue of the operator Wa- However, this operator does not 
conserve probability [that is, sets the time dependence of fU(C,s,t), but the 'total probability' ^^Pyi(C,s,t) is 
not a constant of the motion, except at s = 0]. To interpret this non-conservation, we define a new stochastic process 
(a 'modified dynamics') with rates Ws{C' — > C), chosen so that we can decompose (|16p as 

dtpA{C,s,t) = [WsiC ^ C) - rsiC)Sc,c'] PA{C',s,t) + 6rs{C)PA{C,s,t) (61) 
c 

with rs{C) = Ws{C C) and 5rs{C) ~ ?'s(C) — riC). This decomposition is discussed in appendix [Bl and the 
rates Ws{C — > C) are given in (|Bip . 

For the purposes of the cloning algorithm, we note that the first term in (j61|) conserves probability (in the sense 
given above), while the second term represents the creation or destruction of copies (clones) of the system. That 
is, starting from a large number of copies of the system, we let each copy of the system evolve with the modified 
dynamics (rates Ws)- In addition, the copies are subject to a creation/destruction process with a configuration- 
dependent rate 6rs(C). In this way, the number nc\{C, t) of copies of the system in configuration C at time t has the 
same time evolution as -Pa(C, s, t) in To avoid the ensuing exponential increase or decrease of the total number 

of copies [which behaves as e*'''-*'^*''], one compensates the clone creation/destruction rates of (|6T|) with configuration- 
independent creation/destruction rates. The rates are adapted as the simulation proceeds, in order to keep a constant 
clone population [l^, [23| ■ These adaptively determined rates can then be used to obtain the dynamical free energy 
^a{s). 

2. Results 

Using the cloning method, we investigated two classes of KCMs. In Fig. 2] we consider spin- facilitated models: the 
FA model in one dimension and the East model in three dimensions. We evaluated the free energy density j^ipxis) 
for various system sizes. Its behavior as a function of iV drastically depends on the sign of s, as is also the case for the 
order parameter pk{s). Negative values of s correspond to active histories, with non-zero mean density of particles, 
while for positive values of s, the mean number of particle in the system remains finite, leading to a zero density and 
activity in the infinite system size limit. 

In Fig.O we consider two models with particle conservation: the KA and (2)-TLG models, both in two dimensions. 
Remarkably, the picture is very similar to the previous one, the (conserved) density being replaced with the order 
parameter rx(s). 

In the active phase (s < 0), the order parameters pk{s) and rx(s) converge rapidly as the system size N is increased, 
for all the KCMs that we considered. On the other hand, in the inactive phase (s > 0) the order parameters decrease 
with system size as . Comparing Fig. [Hand Fig. [2] confirms the analysis based on variational bounds on iPk{s)- in 
the limit of large system size, KCMs exhibit dynamical first-order transitions at s = 0. For models on finite lattices, 
the equilibrium (s = 0) dynamics are representative of the active phase, and the system crosses over to the inactive 
phase at a value of s that scales as for large N. 
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FIG. 4: Finite-size scaling in the s-ensemble associated to K, in the vicinity of the singular point Sc — for KCMs with 
non-conserved number of particles: Fredrickson- Andersen model in dimension d — 1 (top) and East model in dimension d = 3 
(bottom). The temperature is T = l/f3 — 0.91, and linear system sizes L are given on the graphs. The finite-size scaling 
illustrates the first order dynamical phase transition in Sc = 0. (Left) large deviation function jj-ipKis). (Right) density of 
excited sites pk{s). 



D. Criticality at zero temperature and dynamical phase transition 

We emphasize that while a zero temperature dynamical critical point is common to many KCMs (36j . this is not 
a sufficient condition for dynamical phase coexistence. Rather, the relevant feature is the presence of states with 
subextensive escape rates, as discussed in Section llV Al In this section, we consider the AA model. The FA and AA 
models both have zero-temperature dynamical critical points, with the same scaling exponents and closely related 
correlation functions 37|]. However, all states in the AA model have extensive escape rates, so we do not expect any 
transition at s = 0. In the following, we show that this is indeed the case, by discussing the AA model both within a 
mean field approximation and using exact results in one dimension. 



1. 'Mean-field' variational hound 



We consider the bosonic AA model in dimension d. Following Section flVBl we calculate the variational Landau 
free energy using the Doi-Pcliti representation, obtaining 

Tk{p, s) = 2d[2Dp{l ~ e-') + fc' + kp^ - ie'^p^/kk'] (62) 

and we identify '0var(s) = — miup s] as a lower bound on Nipxis)- 

The variational estimate for the steady-state density and the variational bound are 

p_(s)=e-^y^+(c-«-l)^ (63) 

V'var(s) = 2d(fcpL(5)-fc') (64) 
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FIG. 5: Finite-size scaling in the s-ensemble associated to K, in the vicinity of the singular point Sc — for KCMs with 
conservation of particle number: Kob- Andersen model (top) and Triangular Lattice Gas model (bottom) Particle density is 
p = 0.5, and linear system sizes L axe given on the graphs. The finite-size scaling illustrates the first order dynamical phase 
transition in Sc — 0. (Left) large deviation function j-j%I)k{s). (Right) order parameter j^(r)a- 
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FIG. 6: Dynamical free energy density ^%1}k{s) for the AA model in dimension one. It presents no singularity at s = 0. 
Although the AA model can be mapped to the FA model and displays the same critical properties at zero temperature, it is 
not subject to dynamical phase coexistence. 



The variational bound is is indeed analytic for all s, consistent with our intuition that the A A model has no dynamical 
phase transition. Again, these variational estimates are exact for the mean-field AA model in the limit of large system 
size, if we set 2d= 1. 



2. AA model in one dimension 



In addition to the mean-field case, we can also obtain the large deviations of the AA model in one dimension, 
through a mapping to a free fermion system. The evolution operator associated to K for AA model can be written 
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in a spin-i representation (recall Eq. ((35)) and see, for example, Ref. [46|): 

i 

- kfiifii+i - k'{l ~ nj)(l - hi+i) - Dhi{l - n^+i) - Dni+i{l - (65) 

where z = e~*, af — ^{af ± icrf), fn — ^-^^ and af,a'^,a^ are the usual Pauli matrices. In the spin language, the 
presence (or absence) of a particle at site i is coded by an up (or down) spin. We use the detailed balance property 
to symmetrise this operator followed by a Jordan- Wigner transformation |46l . l47j 

cr+ = a exp (iTi- ^ cj Cj ) a- = exp («7r ^ cj Cj ) c| (66) 

which allows us to represent the spin operators in terms of fermionic creation/annihilation operators cj and cj. For 
values of the parameters verifying k + k' = 2D [i^ , this puts Wk into a quadratic form: 

W_R- = — ^ ^{k — k' — z{k + k') cos g)cjcg — izy/kk'cl^c^^^ sin q + izVkk'c^qCq sin q — k' N (67) 



where we introduced Fourier-transformed operators Cq = Cje""^^ and = cje~*^^ . 

We now introduce new fermionic operators Pq = cosOqCq — isinOqcLq, (3^ = cos^gcj + is'mOqC^q. Taking 7r/4 < 
9q < 7r/4, we write sin20g = {2z^/kk' sing/Jig), with 



nq = ^y4:kk'{z^ - 1) + [k + k' - z{k - k') cos(7]2 (68) 
Then, the dynamical free energy is the largest eigenvalue of the operator 

Wx = ^5][f^,(l-2/3t/3g)-(fc + fc')] (69) 



2 

Finally, for large N, we convert the sum over q to an integral, arriving at 

N 



(fc + fc')+ / 



(70) 



which depends on s through the dependence of 57, on z = e"**. 

This exact result for the AA model in ci = 1 shows that the large deviation function iPk{s) is analytic at s = (see 
also Fig. [6]), as opposed to the one of the FA model. Despite the presence of a dynamical critical point, the AA model 
has no configurations with subextensive escape rate r{C), and does not exhibit dynamical phase coexistence (in the 
vicinity of s — 0). 

V. PROPERTIES OF TRAJECTORIES IN THE s-ENSEMBLE 

We have proven the existence of a first-order dynamical phase transition in KCMs, and compared the behaviour 
of these models with the A and AA models. The effect of the field s is to generate an ensemble of histories, biased 
towards small or large activity. In order to gain insight into this transition, we now discuss the histories that dominate 
the s-ensemble when s is finite. 

A. Effect of temporal boundary conditions in the s-ensemble 

1. General considerations 

In steady states, the (unbiased) ensemble of histories is invariant under translation in time. Suppose that b = b{C) 
is a configuration-dependent observable, and B = At' b{C{t')). Then, for trajectories of length t the expectation 
value of the observable b at time r is 

{b{T))^\{B), (71) 
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independent of the time r. 

However, introducing a field s biases the ensemble of histories, and, in general, time translation invariance is broken. 
This effect is a dynamical analog of boundary effects in classical thermodynamics: if a system is finite, the behaviour 
near its boundaries is different from that of the bulk. In the s-ensemble, we consider trajectories C(t): the boundaries 
of the trajectory are r = and t — t, while the analogy of the 'bulk' is <C r <C t. In the limit of large time, extensive 
quantities are dominated by the bulk: we have 

(6(r)>,^l(i?)„ (72) 

for <^ T <^ t. However, in general we have (6(t))s ^ {b{t))s 7^ (^(0))s- (In section lV A2l we illustrate these differences 
by calculating (6(r)) in the A-model.) 

More generally, it is possible to express the average at the final time, {b{t))s and the time average 7(5)^ by 
considering the eigenvalues and the eigenvectors of the operator Wa (see appendix . Using this approach, one can 
perform a perturbation theory around s = 0. In particular, when detailed balance is verified and s is conjugate to an 
observable of type B, the bulk and boundary averages differ at first order in s: for large times 

{b{t))s = (5>+,s6« + ... (73) 
{B)s - i [(6) +25^1) + ...] , (74) 

where an explicit expression for fo'^-* is given in Eq. (jA13p . 



2. Effects of temporal boundaries in the A-model 

We now illustrate the effect of temporal boundaries in the s-ensemble, using the (bosonic) A-model. We define the 
average particle density in this ensemble at a time t: that is, 

Pis;T,t) ^ ^{ntot{T))s, (75) 

with < T < t and n^ot = X^j^i- To completely specify the problem, we must set the initial conditions in ([7]): we 
take a Poisson distribution with mean density cq: 



Po 



(Co) = n^ (76) 



To proceed, we write A/Vihist] = j^M-^^t) (j^'j^^^^.^^')^ define P {iitot, J^t , K,t) to be the probability that the 

system contains ritot excitations at time t, having made K changes of configuration, and with the observable A/V [hist] 
taking a value A/V- Then, we define the generating function 

P„,„, = P{ntouh,s,t) =Y. I cW;c-^-^^-^^P(n,A;,if,t) (77) 

K •' 

SO that 

p{s;T,t) - -___ln^P(ntot,/i,s,0 (78) 



Deriving an equation of motion for Pntot is a straightforward generalisation of the derivation of (|16p : the result is 

9tP„.o. =e-^[cArP„,„,_i + (n+l)P„,„,+i] - [cA^ + n + /ine(r-t)]P„,„, (79) 

where 8(i) is the usual Heaviside step function. With the initial condition of (fTG]) . this equation of motion is 
solved by a Poisson distribution with a time-dependent normalisation factor: we take Pijitot, h, s,t) — exp[^'(i) — 
Po(i)]Po(0"'°V(^*totO- Then, the mean, po{t), and normalisation factor, 'i'{t), obey 

poit) = cNe-' - [1 + heir - t)] po{t) (80) 
i'it) = po{t)+e-'po{t)~cN (81) 
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p{s; T, t) 



intermediate regime 




0.1 initial transient regime final transient regime 



f T 

5 10 15 20 25 30 

FIG. 7: Mean density as a function of time in the A-model, for the s-ensemble of histories of length t = 30. We take Cq = 0.1, 
c = 0.2, s = -0.5. 

with initial conditions po(0) = cqN, ^{0) = 0. We identify e*(*) — ^ P{ntot, h, s,t), so we solve for ^'(t) and 
use ([751) to obtain 

p(s; T, t) = cz^ + e"*(l - z)(co - cz) + e~'' z{cq - cz) + e^"*cz(l - z) (82) 
where we defined z = c^'^ for ease of writing. 

The r-dependent density p{s;T,t) exhibits four different regimes (Fig. [71). 

• Short trajectories. For r < t ^ 1, the system has a density p{s; r, t) ~ cq close to the density at time 0. 

• Long trajectories, stationary (bulk) regime. For 1 ^ t <^ t, the system adopts a density p{s; r, t) ~ ce"^**, inde- 
pendent of r. For long trajectories, this average value coincides with the time averaged density dTp{s; r, t). 

• Long trajectories, initial transient regime. For early times r <C 1 ^ i, the density depends on the value of s. 
This dependence persists even for r = 0: that is, the trajectories that dominate the s-ensemble have non-typical 
initial conditions as well as non-typical bulk properties. To be precise, /5(s;0, i) = coe~*: the influence of the 
initial condition decays into the bulk as p{s; t, t) ~ (cqC"* — ce~^^)e~'^ + ce~^*. 

• Long trajectories, final transient regime: for t ^ t the density at the final time t is p{s;t,t) = ce^^. Moving 
away from this boundary, the density decays into the bulk as p{s; r, t) ~ (ce~^ — ce~^^)e'^~* + cc^^**. 

We note that if the initial density cq is equal to the equilibrium density c, then the ensemble at s = has time- 
reversal symmetry. Since the observable K respects this symmetry, the s-ensemble is also time-reversal symmetric 
p{s; r, t) — p(s; t — T,t), for c = cq. In this case, the initial and final transient regimes are related by this symmetry. 

We have used the A-model to calculate the time dependence of p{t) exactly. However, we emphasise that the four 
regimes identified here are very general. When t is large, histories in the s-ensemble are characterised by an extended 
intermediate (bulk) regime, with initial and final transient regimes that decay exponentially into the bulk. 

B. Landau-like theory for fluctuations within the s-ensemble 

In this section, we study large deviations of observables within the s-ensemble. For example, for an s-ensemble 
parameterized by the observable K, we consider the probability of observing a history with a particular value of an 
observable B unrelated to K. In particular, we connect the large deviations of the average excitation density p to the 
variational free energy J-k{p, s), defined in (|55p . 

1. Variational calculation ofipKis) in a general mean-field model 

We consider systems for which we can write the master-like equation (|16p in terms of a single co-ordinate ritot- In 
the mean-field FA model, this co-ordinate is the total number of excitations, but it might also represent (for example). 



19 



the total magnetisation of a mean- field Ising model [8|. To be precise, we assume the master-like operator Wk has 
matrix elements 

(Wk)„,„, = e-^W+_, Sn',n-i + e-^W-^, 5„^„+i - [W+ + M/-]5„-,„. (83) 

where we have abbreviated ritot to n, for compactness, and are the rates for transitions from the state n to state 
nil. For an example, consider the (bosonic) mean-field variant of the FA model, for which = W{n n + 1) = cn 
and W~ = W{n ^ n — 1) ~ (n — l)n/N, as defined in Section Fill CI We have assumed for convenience that all 
processes in the system change the co-ordinate ntot by one: the generalisation to other cases (such as the mean-field 
AA model) is straightforward. 

Following Section [II B 6[ we now symmetrise the operator Wk, so that the dynamical free energy 4'k{s) is given by 
the largest eigenvalue of the operator 

(Wk)„.„' = (W^+M/„-+i)'/'e-^,5„,,„+i + (M^+M^„-+i)i/2e-«<5„,„,+i - iW+ + W-)5^,n' (84) 

For large systems, {N co) the eigenvector associated with the largest eigenvalue takes the form Vn — e"^-''''''/^ 
with p ~ ntot/N, and the function /(p) has a unique global minimum. Then, Eq. (|27p states that 

E„,„,{e-lW^+„,l^„-J^/^(e-/'('''+e+/'(''))-l^„t,-l^-4e 
M^) - max ^^^^^^e^^/(p) ^^^^ 

For any trial function /(p), the sums over ritot in (|85p are dominated by the occupation numbers ntot such that /(p) 
is minimum (which implies in particular /'(p) = 0). Thus, the direct dependence on / vanishes: we are left with a 
maximisation over the position of the minimum in /(p). Since the form of /(p) is irrelevant, the eigenvector Vn can 
be written in the form given in Eq. (j54p . Using this choice, wc arrive at 

4'k{s) = - min (p, s) (86) 
p 

where the variational free energy J-k{p,s) was originally defined in (|55p . For these mean-field models, it takes the 
form 

As discussed in Section llVB[ J-k{p, s) gives a bound on tpKis) for all systems. However, for systems with mean- field 
geometry, we showed that the form of the trial distribution is irrelevant in the limit of large system size. Thus, we 
write Eq. (|86p with an equality, and not as a bound. We now discuss the physical interpretation of this result. 



2. Physical interpretation of the variational free energy 

Consider an s-ensemble in which trajectories are weighted by the usual factor e^^^I^'^'l, but with the further 
restriction that the time-integrated density be fixed. That is, we write the (unnormalized) probability, in the s- 
ensemble, to measure a time-averaged density p. 



where the asymptotic behaviour of the left hand side at large t defines the function J-^{p, s) 
Taking a Laplace transform of ((88|) with respect to p, we arrive at 

{y'i^{p,s)+hNtp) 



Zkm{s, /i) = ^ exp - - hNj )= j e''^^'<^P''>+'''''P> (89) 

where we write N = JJ" dr ntot {t) , noting the similarities with the generating function of (j78p . 

Now, by analogy with (jl4p . we identify Zx j^{s,h) as the partition function for an '(s, ft.)-ensemble', in which 
histories are biased both by their activity K and their time-integrated number of excitations Af. Repeating the 
analysis of section Hi B 3[ we observe that 

'ipK,j\r{s,h) ^ lim \ \nZK,Ar{s,h) (90) 
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is the largest eigenvalue of an operator W/^.aa, whose elements are 

(WK,AA)n',„ - W;,ti e-^5„',„-i + e~'6n',n+i - [W+ + W- + hn] ^„,,„ (91) 

The largest eigenvalue of this operator can be obtained by symmetrising and repeating the variational analysis of 
the previous section. The result is 

V'-ft:,AA(s,^) = -niin(J^if(p, s) + /ip) (92) 

p 

which applies in the limit of large system size N [since in that case, the maximisation over the function /(p) can be 
replaced by a maximisation over the density p]. However, performing a saddle point analysis directly on (|88|) reveals 
(for large times t and finite system size N), 

fpK,Af{s, h) = - min(:r^(p, s) + hp) (93) 
p 

Thus, in the limit of large system size N, the Legendre transforms of J-k{p,s) and J-'^{p,s) are equal. It follows 
that the large deviation function T'^{p,s) coincides with the variational free energy Tk[p,s) as long as the inverse 
Legendre transform can be performed, However, in Section [IVBI we showed that in KCMs, J^{p, s) typically has two 
minima, separated by a range of densities in which it is 'non-convex ': dlT{p,s) > 0. In this case, the inverse Legendre 
transformation cannot be performed. In fact, the non-convexity of J- {p, s) arises because histories with some values 
of p are unstable in the s-ensemble, as we now show. 



3. Non-convex free energies: phase separation in time 

In the thermodynamics of finite-dimensional systems, one typically has s" (e) < o[2l|. Loosely, this property arises 
because any energy density e can be achieved by separating a system into two regions, separated by an interface whose 
energy cost scales subextensively with the size of the system. Thus, the total energy density is e = (1 — x)ei + xe2 + 5 
where ei and 62 are the energy densities of the two regions, x is the fraction of the volume of the system taken 
up by the second region, and 5 is the energy of the interface divided by the total volume, which vanishes in the 
thermodynamic limit. This leads to the 'lever rule' e = (1 — x)ei + xe2- The total entropy density associated with 
these configurations is s(e) = (1 — x)s{ei) -\- xs{e2), and using the lever rule, it follows that s"(e) < 0. However, in 
mean-field geometries, interfaces cannot be formed, and this argument cannot be applied. 

Interestingly, in the statistics of histories, phase separation is possible even in mean-field systems. We consider 
the large deviation function J-^{p, s), at a density p for which J-k{p, s) is non-convex. We will find that the average 
in (|88p is dominated by histories that phase separate in time. To prove this, we use the methods of Donsker-Varadhan 
theory, described in appendix [C] This method allows us to prove that 



IV) {V\V) 



(94) 

{V\p\V)=p 



where the minimisation is over distributions V{C) such that ¥{0)^ p{C) — p. By analogy with the thermodynamic 
case, we take V{C) = (1 — x)Vp^ (C) -I- xVp^ (C), where Vp{C) was defined in We then minimise over the densities 
pi and p2, choosing x = {p — pi)/{p2 ~ Pi) to ensure that the mean density is p. Taking x = 0, we have a bound 
T^{p,s) < !Fk(Pts). However, if J^k{p,s) is non-convex (that is, dlTK{p,s) < 0) we can find a lower bound on 
J^^(p,s) that is smaller than !Fk{p,s). For example, in the FA model in finite dimension (and in the limit of large 
system size N), we find that J-'*{p,s) is minimised by pi = 0, p2 = cz^ and x — /3//O2, for < p < ce~^*. For 
< p < cz^ ^ this variational approximation to J-^{p,s) indicates that the system separates into two phases with 
densities and cz^. We arrive at a bound 

^KiP^s) ^ f pc(l-z^), p<cz^ , . 

-1 n(^J-n-r>^ fTTTW ^^^^2 lyOj 



2d ~ \ P(c + p- 2z^/cp), p>cz^ 

from which we note that dp!F*{p,s) = in the two-phase regime: this is the Maxwell construction [2l[. The result 
for the mean- field FA model is obtained by setting 2d = 1, as in Section flV Bl in that case, the bound is saturated 
(this follows since J-* is convex and its Legendre transformation is known to be equal to that of JF) . 

In addition to establishing the convexity of (p, s) , Donsker-Varadhan (DV) theory also provides an interpretation 
of the distribution V*{C) that minimises (|M)) . (We normalise to J2c^*i^)^ ~ ^ ^'-^^ convenience.) At large t, we 
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consider histories in the s-ensemble with fixed average density p. The DV theorem states that this sub-ensemble 
is dominated by trajectories for which the fraction of time spent in configuration C is n*{C) = V*{C)'^. If the 
distribution ^*{C) is dominated by configurations C with density p, we conclude that the histories in this sub-ensemble 
are homogeneous in time. However, if fJ-*(C) is associated with a bimodal density distribution, it describes histories 
comprising separate periods of of time, some with low excitation density pi and some with high excitation density p2- 
Finally, we give the interpretation of Tk{p, s). In systems described by the single co-ordinate ntot, the DV theorem 
states that ^-^^^^{9,5) jg ^-j-^g (unnormalised) probability (in the s-ensemble) of a history in which almost all config- 
urations have density ijitot/N) equal to p. This can be compared with the probability e~'^*--'^'<ip^'^) of a history with 
a time-averaged density p. In this sense, Tk{.P^ s) can be interpreted as a Landau- like free energy for homogeneous 
trajectories, while T^(p,s) is the large deviation function for the density p. In finite dimension, Tk{Pts) provides a 
bound on the large deviation function !F'^{p, s), based on the assumption that histories are spatially and temporally 
homogeneous. That is, fluctuations in space and time are neglected. In Section FV C[ we discuss how these fluctuations 
can be taken into account. 



4- Landau-like free energy in other s-ensembles 



So far, we have considered the large deviations of the density p in an s-ensemble that is defined in terms of the 
activity K. The variational free energy can be simply extended to s-ensembles defined as in (I12p . Consider again a 
mean- field model specified by rates and an observable A of the form given in Q , which is incremented by for 
transitions from state n to nzL 1. Repeating the analysis of Section IVB H we find iPa{s) ~ miiip !Fa(p, s) with 



^a{p,s) = -{-2iW+^e 



(96) 



For example, in the case of the complexity Q^, one has — In ■ 



— and 







1-s 




tp+{s) — — min < —2 
P { 




2 


W^p + W^^p 



W^p + Wj,^ 



(97) 



where, again, this variational bound is exact because of its independence of the form of the trivial wavefunction used. 
By analogy with Tx{PtS), we find that e~^^^A{p,s) gyve the probability of homogeneous histories with density p, in 
the s-ensemble. 



C. Dynamical free energy landscape (field theoretic approach) 

In Section [VB| we have considered large deviations of time-averaged observables, using these quantities to charac- 
terise the histories within the s-ensemble. We now discuss the calculation of dynamical correlation functions within 
this ensemble. We make use of a field-theoretic description of the FA model. 



1. Field-theory for the bosonic FA model 

Using the Doi-Peliti representation of the bosonic FA model (Section IIII C|) . we use coherent states to write the 
partition function Zif(s,i) as a path integral over (time-dependent) functions {<pi} and {ipi} [11]. Then, taking the 
continuum limit, we promote these functions to fields {4>xTT(t>xt) depending on position x and time r, where (j}xT has 
the dimensions of a density and i/)^,- is dimensionless. Introducing sources h and h for the fields 4> and 4>, we write 



Z[s, t; hxr, Kr] = J T^[4>xt,4>xt] exp | - S'/f [0, + y dxdr {hxr<PxT + hxAxr) \ (98) 

where the path integral is over histories of duration t, and (see, for example [sj) 



SK[4'xr,4'xT] = / dxdr 



\^^xrdt(t>xr - Ml^ {$xr<l>xr+do'^)-e~'{(j)xr+clQ'^$xr) (1 + V^)<^:.r 0xr } (99) 



where /q is the lattice spacing, and we have taken a gradient expansion, truncating at quadratic order. We identify 
ZKis,t) = Z[s,t;0,0]. 
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2. Saddle point approximation 

We now show that a saddle-point analysis on the action recovers the results of the previous sections. The saddle- 
point equations are obtained by maximising the action with respect to (j) and 4> in the absence of the sources (h, h). 
The saddle occurs for fields that are homogeneous in space and time, with magnitudes satisfying 

= 24>(j){(t)l'},-e-''c)+(t){c-e-'(t)lf^) , (100) 
= 2 00 + C(^ (l - e""^) . (101) 

These reduce to a single equation if we take cIq^(}> = (f> (this origin of this symmetry becomes clear if we use the 
symmetrised operator W in the construction of the original path integral). In this single variable, the solutions are 
6 = and 



+ lV9p^^^ (102) 



(p = (pact = . ^ 

To estimate the dynamical free energy, we simply identify iPk{s) with (— t"^ minS'[(/), (j)]) where the minimum is over 
value of the action at the two saddles. The result is 

^-(^)"{L(c<^.et)^[0acte--l], T<'0 (^03) 

where the approximate equality indicates that we are working in the saddle-point approximation. We identify the 
time-dependent density (per site) of excitations in the s-ensemble {p{t))s = ^o(0*(r)(/)(T)). Away from temporal 
boundaries, we take the saddle point value for this average, obtaining 



It is easily verified that (|103|) and (jl04p coincide with the variational estimates ([58|) and (|59p . 

In principle, we can can now use the tools of dynamical field theory (49j to incorporate fluctuations around the 
saddle points, and to calculate spatiotemporal correlation functions in the s-ensemble. For example, defining a density 
field n(x,T) through a contimmm limit of the original occupation variables n^, we have 



{n{x , T)n{y , t' )) s = liiZ[s,t;h,h] 

6h{x, T)Sh{y, T')Sh{x, T)6h{y, r') 



(105) 

h=h=0 



Thus, for models with a field-theoretic representation (such as the FA model), the framework described in this 
section provides methods for systematic calculation of correlation functions and fluctuation effects in the s-ensemble. 
However, these field-theoretical calculations beyond the scope of this paper. We emphasize that the analysis of 
Sections IIV Al and IIV CI establishes that a dynamical first-order transition does occur at s = in finite-dimensional 
KCMs. Thus, while we expect fluctuations to have quantitative effects, the qualitative picture obtained through this 
saddle point analysis is not changed. 



VI. OUTLOOK 



We have analyzed the dynamics of kinetically constrained models, using an ensemble of histories which span a long 
time t. This analysis used dynamical tools 0, 0, 0, H, [lH constructed by analogy with the usual Boltzmann-Gibbs 
theory of equilibrium systems. We have established that this procedure captures physically relevant features that are 
not accessible from the steady state distribution of configurations in these models. 

We have shown that the steady state of KCMs lies on a first-order dynamical transition line, characterised by a 
coexistence between active and inactive histories. This first-order line is present both in mean-field systems and in 
finite-dimensional models. Its existence is proven by variational bounds on the dynamical free energy, and confirmed in 
numerical simulations of several kinetically constrained models, including both spin-facilitated models and kinetically 
constrained lattice gases. We have defined dynamical Landau-like free energy, whose form is intimately connected to 
the existence of dynamical heterogeneities. 
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Earlier studies of non-equilibrium systems used a similar thermodynamic formalism for dynamics to reveal first-order 
transitions arising from a static phase transition [50l[5l| or from an absorbing state [8]. To place our work in context, 
we emphasise that our dynamical phase coexistence is not related to such phenomena. However, the transitions in 
these models all appear as singularities in their large deviation functions, consistent with the idea |52l | that phase 
transitions both in and out of equilibrium can be studied through the eigenvalue spectra of their master operators. 
Moreover, the focus of the current paper is on transitions between stationary, time-reversible dynamical states, and 
therefore we concentrated on large deviations of quantities that are symmetric in time: this is to be contrasted with 
studies that have concentrated on currents of entropy or particles j42, j53] , although recent work has hinted that large 
deviations of time-symmetric observables may also be of importance in non-equilibrium steady states [s^ l . 

We expect our approach to be meaningful in a wider class of systems than those probed in this paper. For example, 
glass-forming liquids are known to be dynamically heterogeneous, and this feature can be captured in computational 
simulations of atomistic models. It would be interesting to establish whether this heterogeneity is linked to a dynamical 
phase transition similar to that present in KCMs. This could indicate a more general link between glassy properties 
(not necessarily related to dynamical heterogeneity) and dynamical phase transitions. 

Finally, we observe that an experimental scheme for sampling the s-ensemble would be very valuable, since it would 
provide a direct test for the existence of a dynamical phase transition. However, the fact that the generalised master 
operator Wa does not conserve probability makes the search for such a scheme rather challenging. 



Acknowledgments 

VL would like to thank Julien Tailleur for useful and continuous discussions, and RLJ and JPG thank David 
Chandler for extensive discussions. This work (EP and FvW) was supported by the French Ministry of Education 
through Grant No ANR-05-JCJC-44482. JPG was supported by EPSRC under Grant No. GR/S54074/01. While at 
Berkeley, RLJ was supported initially by NSF grant CHE-0543158 and later by Office of Naval Research Grant No. 
N00014-07-1-068. VL was supported in part by by the Swiss ENS, under MaNEP and division H. 



APPENDIX A: AVERAGES IN THE s-ENSEMBLE AND EIGENVECTORS OF Wa 



1. Eigenvectors of Wa 

In this appendix, we discuss some properties of the operator Wa, and their consequences for averages in the s- 
ensemble. We write in terms of its left and right eigenvectors |L„) and |i?„): Wa = J2n ^n\Rn) {Ln\ with 
eigenvalues Aq > Ai > • • •■ The maximal eigenvalue Aq is equal to ipAis). One can normalize eigenvectors so that 

(i„|i?„) = (5„„ and (-|i?o) = l (Al) 

where | = J2c(^\ ^^'^ projection state. 

Thus, for long times, we have e'*'-** — |i?o)(-^o|e*''''*^*'' -I- . . . where the omitted terms on the right hand side are 
exponentially smaller than the dominant first term. Therefore, starting from an initial state \Po) = X^c Po{C)\C), with 
(^l^^o) = 1, one has, for large times 

\P{t)) = e'^-'lPo) ~ \Ro)e^-^'^\Lo\Po) + ■■■ (A2) 

where we write the largest eigenvalue of Wa as Aq = ipAis), and the omitted terms are exponentially smaller than 
the first one, for large times t. This allows us to identify the largest eigenvalue of Wa with the dynamical free energy 
limt^aa log ZA{s,t), through Equ. 



2. Time averages 

We now consider a configuration-dependent observable b{C), and an s-ensemble defined as in (jl2[) . using an observ- 
able A of the form given in We provide a link between the eigenvectors of Wa and two weighted averages: the 
average of b at the final time t in the s-ensemble 

{Kt))s = j—ji (A3) 
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and the time-integrated average of b in the s-ensemble 



-sA 



f^dTh{C{T))c 

iB)s ^ ^ '- (A4) 

As discussed in Section fV Al {B)s grows Hnearly in time, but, in general dt{B)s ^ {b{t))s. In operator notation, 
the average of b at the final time is 

,,,,,, (-|&e*^-|Po) 



where b denotes the diagonal operator of elements 6(C). Using the normalization (IAl|) . and the large time result (jA2p . 
we arrive at 

= Hb\Ro) (A6) 

Thus the right eigenvector |i?o) gives the distribution over configurations C at the final time t. 

On the other hand, the integrated average (i?)s is obtained from the mean value (6(t))s in the intermediate regime 
< T < 

For ^ T <C i, we have 

e^^-|Fo) = e^'^-(^)|i?o) (io|Po) (A8) 
(_|e(*— )Wa ^ ^_|p^j^ (Lo|e(*-")'''^(") (A9) 



and hence 



Ub)s = {b{T))s = {Lo\b\Ro) (AIO) 



t 

Thus, while the average {b{t))s depends only on |i?o), the average {B)s depends on both |i?o) and {Lq 



3. Dynamics with detailed balance 

From it follows that if a system obeys detailed balance, its master operator satisfies = P(,~^WxPoq, where 
Poq is a diagonal operator with elements Poq(C). Thus, = Pp~^|P„). Using this property together with results 
from the previous section, and denoting |Po) = J2c ^o(C, s)\C) we write 

= ^6(C)Po(C,s) (All) 

c 

d,{B). = T.biC)^fj^ (A12) 

Clearly, these averages are not the same in general. Expanding about s = and using Rq{C,0) ~ Pcq(C), we arrive 
at dini), with 

b(^)=ybiC) ^^f''^ (A13) 
c 

Finally, we note that expectation values of the form dt{B)s take a simple form when written in terms of the 
eigenvectors \V) of the symmetric operator , discussed in Section Hi B 61 The matrix elements of this operator are 

(Wif)c,c' = Pcq'^'(C)(WK)c,c'PcV'(C), so its eigenvectors are F„(C,s) ex V^„(C, s)P„(C, s) = Rn{C , s) / ./iQfi) . 
Thus, we have (for large time) 

hs). ^ ^^^f (A14) 

which links the eigenvector Vo{C,s) to physical observables such as B. 
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APPENDIX B: OBSERVABLES OF TYPES A AND B 



Here, we discuss the connections between observables of the forms given in ([9]) and ([22]) : we refer to these observables 
as type A and type B respectively. We begin with a result that is used in the numerical methods of [H, [lO] ■ 

Consider an s-ensemble defined as in Section fll B 31 That is, take a system with rates W{C — > C) and modify the 
statistical weights of its histories by a factor e~^"^, where A is an observable of type A. In addition, we define a second 
stochastic process ('modified dynamics') through the transition rates 

W^,(C^C') = e-^"(c,c')yy(c^C') (Bl) 

where the a{C,C') are obtained from the definition of the observable A, through In addition we define two 
configuration-dependent observables, rs{C) = J^c ^s{C C), and 

Srs{C)=rs{C)-r{C) (B2) 

[Here, r{C) — J^c C) is the escape rate for the dynamics W , as in the main text.] 

Motivated by the decomposition of Equ. we can establish two ways of defining the same s-ensemble. From (O, 
it is easily verified that 

Prob[hist|W]e-*'^ Prob[hist|IFs]e''^= (B3) 

where the notation Prob[hist|Vl^] refers to the (unmodified) probability of a history in a system with dynamical rates 
W, and 

5Rs = I dT(5r,(C(T)) (B4) 

Thus, Equ. (jBSp states that histories in the s-ensemble parameterized by A for the original dynamics W have the 
same weight as histories in an s-ensemble parameterized by (5i?s, for the modified dynamics Ws- 

Since the two ensembles are identical, it follows that all observables have the same averages: for example 

{Oe-'^)w ^ {Oe'''^)w., (B5) 

where the subscript on the average refers to the dynamical rules used for the sum over histories. Further, this result 
holds for histories of finite duration as long as the same initial conditions are used in both averages. 

For the specific case where the observable A is the activity K then this relation takes a particularly simple form. 
Following Section fll B 31 with a(C,C') ~ 1 for all C and C, the master operator associated with this s-ensemble has 
matrix elements (W^)^ ^, = e"" W{C' -> C) - r{C)5cfi' From (|BT|, we find Ws{C C) = e-'W{C C): that is, 
the modification to the dynamics simply involves a rescaling of time by a factor e~*. In addition, for B = R, we have 
Srs{C) = sr{C), so we define an s-ensemble associated with the observable i?[hist] — dr r(C(r)), which is of type B. 
From the analysis of Section fll B 5| the master operator associated with this ensemble, Wr, has matrix elements 

0^r)c,c' = WiC -> C) - (1 -h s)riC)6c,c' (B6) 

from which we can see that 

WKis^^c-'WR^c' -I) (B7) 

This equation relates the dynamical free energies of the s-ensembles for K and R, based on the same unbiased 
dynamics W. The dynamical free energies iPk{s) and (/)i{(s) are given by the largest eigenvalues of and W_r: they 
satisfy 

^if(s) =e-^0fl(e^-l) (B8) 
Hence, we can also relate the cumulants of the observables K and R. For example, 

(K) = (R) (B9) 
mc = {R')c {R) (BIO) 
This last equation provides an interpretation of the variance (second cumulant) of K, through 

-{K'), = -{r) + -[{r{t)r{t')) (r(i)) (r(t'))] dt' (BU) 
where the correlation function is evaluated at s = 0. 
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APPENDIX C: LINK TO DONSKER-VARADHAN THEORY 



As in the main text, we consider a Markov process described by transition rates W{C C) between configurations 
{C}. For a history C(t), we define the experimental measure 

fliC,t)= f drSccir) (CI) 
Jo 

This history-dependent observable simply counts how much time was spent in configuration C between and t. This 
is the central object of Donsker-Varadhan [sHi theory (see also [1, [sHi)- For large times, the experimental measure 
approaches the steady state distribution, limj^oo j{P-{C,t)) = Pst(C). 

Donsker-Varadhan theory gives information on the large deviations of the experimental measure p,{C,t) from the 
steady-state distribution, in the long time limit. Therefore, it is naturally connected to the statistics of histories and 
to the dynamical ensemble approach discussed in this article. For example, consider an observable b{C) depending 
on the configuration of the system. The experimental measure p,{C,t) determines the time-integrated value of the 
observable b through 



f dTh{c{T)) = Y,b(c)ri{c,t)^t{b)-, 

Jo r 



(C2) 



which defines (fe)^: the average of b with respect to the experimental measure Jl. In this appendix, we establish links 
between the s-ensemble approach and the results of Donsker and Varadhan. In particular, we develop a variational 
method that gives the large deviations of an observable S, in the s-ensemble defined for an (unrelated) observable A. 

1. Donsker-Varadhan large deviation function 

The Donsker-Varadhan (DV) theorem [5^] states that in the long time limit 

Prob [/2(C, t) = tfi{C)] = e* -^[^1 (C3) 

with (see for instance [sl-ls^) 

= {^(^ ^ ^'^^ ^(^^ " r{CMC)6c'c} (C4) 

where the infimum has to be taken over normalized measures p(C), with p{C) = 1. 

If W obeys detailed balance with respect to an equilibrium distribution Pcq(C), the infimum is obtained for p{C) = 
^/ p{C) / Pcq(C) and the large deviation function reduces to 

■/oqM =E{[^(C^C')mC'-C)]'/'[Ai(C)A*(C')]^/'-r(C)M(C)5c'c} (C5) 

C,C' 

for normalised measures fi{C) ~ 1. Writing 



we identify 



as the function to be maximised in (l27l). for the case s = 
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2. Dynamical Landau free energy at s = 

We now apply the DV theorem to the large deviations of an observable b{C). Integrating over a time t, we define 
the history-dependent quantity 

Bit)^ f dTb{C{T)), (C8) 
Jo 

As discussed in Section Pi B 4[ one expects the probability distribution of B{t) to behave as 

f]dy„(5 = - e*'"''') (C9) 

for large times t. 

The large-deviation function 7r(6) can be obtained through the Donsker-Varadhan fimctional using B{t) — 
Y,cKC)ii{C,t), so that 

Pi-oh[B{t) = tb] = {S{B{t) - th)) 

dn {S{B{t) - tb) S{t-^fl{C, t) - ^i{C))) 

diJi S{J2biCUC) - b) {6{fliC, t) - tii{C))) (CIO) 
c 

d^i5{{b)^ - b)e"^^^^ (Cll) 

where the average (6)^ was defined in (jC2[) . Here, we have replaced an average over histories (•) with an integral over 
possible realisations of the experimental measure weighted by their probabilities (which are known from the DV 
theorem). In the limit of large time, we maximise the argument of the exponential, subject to a constraint imposed 
by the (5-function. Hence, 

7r(6) = sup J[/i] (C12) 

/i. with 

which for systems obeying detailed balance can again be expressed in terms of the operator W, using (|C6p . 



3. Dynamical Landau free energy for any s 

We now generalise this analysis to the s-ensemble. We note that the values of 'type A' observables [those of the 
form given in ([5])] cannot be obtained from the experimental measure /2(C). To connect these observables to the DV 
approach, we use the results of appendix IB] 

The large deviations of the observable B in the s-ensemble specified by A are determined by 

^lA{s,b)^{6{B-bt)e"'^)w, (C13) 

where as in appendix iB] the label on the average indicates the dynamical rules used to generate the ensemble of 
histories. From (jBSp . we can write 

nA[s,b) ^ {S{B - bt)e-'^)w - {S{B -bt)e^^-)w^, (C14) 

with an observable SRg and rates Ws{C C) given in Equ. (|Bip and (|B4p . 
Now, following the analysis of the previous section, we have 

^Ais^b) = (^c^'^^5{B{t)-tB))^ = J d^Ji5{{b)^ - B) Qt{J[^^\^A+{Sr.),) ^^15) 

where ((5r)p — X^c J[/i|VKs] is the Donsker-Varadhan functional for stochastic process with rates Wg- 

Again, the integral over /i can be evaluated by maximising the argument of the exponential subject to the constraint 
on {b)^, leading to Q.A{s,b) ^ expi7rA(s,fc) — with 

■T^A{b,s) = - sup Ja[a*,s] (C16) 

/J, with 
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with 

JA[/i, s]=J2 {e-*["('^''^')+"('^'''^)l [WiC ^ C')WiC' ^ C)]'/'[/.(C)m(C')]'/' - riC)^iiC)6c'c} (C17) 

C,C' 

where we emphasise that the rates W{C C) arc those of the original (unmodified) dynamics. From (|C6p . we identify 

l™ ,C18, 

as the quantity to be maximised in (|27p for A = K. Moreover, for h being the occupation number n and A the activity 
K, one recognizes in (|C18p the result ([M]) . 

We observe that these results have been derived for dynamics which obey detailed balance, but they are not 
restricted to that situation. For instance, (jC16p holds in general, with 



J^[m,5] = inf g{e-^"(C''^')W^(C->C')^A^(C)-r(C)/i(C)<5c-c} (C19) 



Finally, we note that these results amount to a generalization of the Donsker-Varadhan theorem (|C3p in the s-ensemble: 
for large times, 

{Je-"-^^(^(C,i)-tM(C))^ ^e*-^-^!^'^! (C20) 
with Jyi[Mi ^] given in general by (|C19p . which reduces to (jC17P if Wa can be symmetrised. 
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